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ABSTRACT: Some of the fundamental neutron star parameters, such as 

the mass and the magnetic field strength, have been 
experimentally determined in accreting neutron star 
systems. We review here some of the relevant data 
and the models used to derive useful information from 
them, concentrating mainly on X-ray pulsars. We then 
discuss the latest advances in our understanding of the 
radiation mechanisms and the transfer in the strongly 
magnetized polar cap regions. 

I. INTRODUCTION 

A significant fraction of the known celestial X-ray point sources 
have been recognized as neutron stars in close binary orbit around 
a normal star companion, while in a number of Y~i^ay burst events 
there are strong indications pointing to a neutron star origin, 
binary membership being as yet unclarified. With a few exceptions, 
we shall be dealing with those systems where the neutron star 
appears as an X- or Y~ray pulsar powered by mass infall. X-ray 
bursts, some of which are thought to be thermonuclear flashes on 
low-magnetic field neutron stars, are dealt with elsewhere in this 
volume (1). Unlike the rotation-powered radio pulsars, some of 
which also emit X-rays and yrays, the accretion-powered neutron 
stars, whether pulsing or not, emit their energy almost exclusively 
above the keV range. Among these systems, those that do show 
regular pulsations are the subject of the most vigorous inquiries, 
since they provide a particularly useful laboratory for testing our 
knowledge of nuclear and particle physics, as well as of general 
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relativistic stellar structure, binary evolution, strong magnetic 
field plasma physics and radiative transfer. Experimental 
determinations of the neutron star masses, magnetic field 
strengths, and indirectly radii, have become available, as well as 
spin periods and period decrements, orbital periods, pulse shapes 
and .electromagnetic spectra above a few tenths of TceV, all of which 
provide a powerful tool, as well as challenge, for theories and 
models . 

While the standard interpretation of binary X-ray sources, wherein 
the gas lost from the visible companion is captured and heated to 
high temperatures by the neutron star (2,3,4), has, on the whole, 
proved remarkably satisfactory, there are a number of tantalizing 
ambiguities about most of the processes involved. In its simplest 
version, the model envisages the mass exchange occurring either by 
slow leakage (Roche lobe overflow) or by capture from a stellar 
wind from the companion. In the first case, the matter's angular 
momentum causes an accretion disk to form, or in the second case, 
its random value leads to a quasispherical inflow. This spiralling 
or quasiradial inflow is assumed to be stopped by the pressure of 
the neutron star's magnetic field at the Alfven surface, unless if 
the field has decayed. The accreted matter is threaded by the 
field and is conveyed to the polar caps at ~ free fall velocity, 
where abrupt deceleration and heating occurs either by particle 
collisions, collective effects or radiation pressure. Pulsed 
emission is due to the spin of the neutron star and the anisotropic 
emitting pattern, or periodic occultation by matter at the Alfven 
shell. We list below some of the relevant X-ray data and their 
interpretation, including some of the yray burst information where 
it relates to neutron stars. 

II. PULSE PERIODS AND MASS EXCHANGE MECHANISMS 

The pulsation periods of upwards of 10 X-ray pulsars are known, 
typically to 6 figure accuracy, ranging from about 1 sec to about 
1Q3 sec (5). A strong clue for distinguishing between the two 
possible modes of mass exchange, has been found for those X-ray 
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pulsars where besides the period P also the period decrement P is 
known. Calculations of the torque exerted on the neutron star by 
an accretion disk interacting with the dipole field of the neutron 
star at the Alfven radius (where B^/8ir ^ pv^) give a well defined 

^ T , • 

theoretical relationship between the period change -P and the 
quantity (PL^/^), l being the luminosity (6,7). The observations 
available for 10 pulsars give a fairly good agreement with this 
relationship utilizing a magnetic field typical of neutron stars 
(see Figure 1). A comparable theoretical relationship for the wind 
transfer case, on the other hand, is not readily available, 
estimates (7) suggesting the absence of a fit. This is partly due 
to insufficient development of the wind theory in an aspherical 
potential, where the atmosphere may also be disturbed by X-rays 
from the companion. The interaction of the wind, and alas also of 
the disk, with the magnetic field is made difficult due to 
uncertainties concerning conductivity, turbulence and viscosity, 
e.g., (8, 9) for reviews. On the face of it, it would seem that 
binary X-ray pulsars tend to appear in systems with accretion 
disks. On the other hand, there is evidence for wind transfer in a 
number of X-ray binaries (10) including some X-ray pulsars (e.g., 
4U0900-40 = Vela X-1 ) . This seems consistent with the fact that 
this object falls outside the "disk" relationship of Figure 1. In 
another case (Cen X-3), there is evidence for' a wind, but it has 
been argued that a disk would nonetheless form (11), although this 
is not the case for most wind sources. Other neutron stars may 
exist in disk systems where pulsations are not detected because of 
the beam not intersecting the observer, or in wind cases, because 
of scattering on the surrounding flow. Interestingly, in the y-ray 
burst of March 5, 1979, where the energy source is unlikely to be 
accretion but it is also believed to be a neutron star (12), an 8- 
second period has been detected (13). 

III. PULSE SHAPES AND PULSATIONS MECHANISMS 

A large body of information exists about the pulse shapes of X-ray 
pulsars (see Figure 2) which contain information on the nature of 
the mechanism for producing pulsations. The pulse shapes range 
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from being strongly dependent on energy (e.g.. Her X-1) to showing 
very little change with energy (e..g., Cen X-3). The multiplicity 
of the pulse structure is often also a function of energy, and the 
duty cycles are invariably large (~ .5). Both symmetric, and 
highly assymmetric pulse shapes appear in different objects. There 
ax-e two principal mechanisms for producing pulses, one of them 
relying on the intrinsic beaming of the radiation pattern, caused 
by anisotropies in the radiative transfer or the accretion flow, 
and the other relying on periodic occultation of the emitting 
surface by an opaque screen near the Alfven surface, revolving at 
the pulse frequency. In the first one, an energy dependent pulse 
profile is expected, the cross sections in the strong field of the 
polar caps being frequency and angle dependent, while in the 
second, only in the photoa’bsorption range 1-6 keV is a dependence 
expected, the opacity at the Alfven surface above 6 keV being 
essentially Thomson. In the case of Her X-1, the observation of a 
soft X-ray flux (14) in antiphase with the hard X-ray pulses can be 
readily interpreted as being due to the hard X-rays being absorbed 
in the screen (15,16,17) followed by reemission at lower 
frequencies. Photo absorption in the range 1 < hv < 6 keV would 
explain the sharp pulse edges, but would also predict stronger 
modulation in this range than at higher energies, which is not 
observed (17). On the other hand, recently some detailed 
calculations of intrinsic beam patterns from homogeneous magnetized 
atmospheres have been made (18,19), which indicate that the duty 
cycle, energy dependence and multiplicity of pulses can be 
reproduced, but difficulties remain in explaining the steep edges 
and the asymmetry of some pulses. The observation of a hardening 
of the spectrum towards the middle of the pulse (20), and the 
apparent trend shown by e.g.. Her X-1, A0535+26 and 4U0900-40 of 
the simple pulse structure breaking up into multiple pulses as 
increasingly lower energies are sampled has also been found in 
intrinsic beam calculations (18). Neither the intrinsic beam nor 
the opaque shell models are self-consistent as yet. In the first, 
the atmosphere structure (gradients) remains to be elucidated, in 
the second, the dynamics and the geometry of the opaque shell is 
uncertain. A very valuable observation would be to determine the 



5 

polarization of the pulses as function of phase, since the 
predicted intrinsic polarization (18) is much greater than for the 
shell case (17) at low energies. 

IV. SPECTRUM 

There is also much information on the spectra of X-ray pulsars, 
most objects having a rather flat power law energy index between a 
few keV to several tens of keV, and dropping abruptly above that 
(e.g., 21,22,23). Line features ascribed to Fe are sometimes 
present in the 6-8 keV region (24,25,26). Other lines, ascribed to 
the cyclotron process, will be discussed in Section VI. The Fe 
lines, having widths of a few keV, may originate near the Alfven 
surface (27,28), since they are rather narrow. The continuum, on 
the other hand, must originate near the surface, to explain the 
large luminosities, often near the Eddington value. As an example, 
the spectrum of Her X-1 is shown in Figure 3. While naively the 
power law photon spectra resemble somewhat an optically thin 
bremsstrahlung, the actual situation, even for a homogeneous 
atmosphere is more complex. The polarization structure of the 
magnetized medium introduces a subtle interplay whereby photons are 
created in one polarization and escape in the other (29,30), the 
thermalization length being a function of frequency. A simple 
homogeneous atmosphere, if only bremsstrahlung and coherent 
scattering is assumed, will look like optically thin bremsstrahlung 
only when very low optical depths are assumed, vJhich do not produce 
the required luminosity for a galactic X-ray source. Increasing 
the optical depth increases the luminosity, but the two-polari- 
zation coherently scattered spectrum then turns over at low 
frequencies too early. The inclusion of comptonization (incoherent 
scattering) may be the answer, although existing Monte Carlo 
calculations (31,32,33) with this effect included still turn over 
too early. Analytic and numeric integration of the incoherently 
scattering transfer equations (34,35) suggests, however, that one 
may reproduce a power law down to low frequencies. 
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V. MASS DETERMINATIONS 

This is a crucial parameter from the evolutionary and microphysical 
point of view, which we review here only briefly. A number of 
neutron star masses have been determined using the dynamical laws 
governing the orbital motion in the binary system. The simplest 
case is that when one has i) the doppler delay curve for the 
pulses from the neutron star, ii) the doppler delay curve for the 
companion (e.g., its absorption line frequencies), iii) an eclipse 
of the neutron star by the companion. i) yields ajj sin i, the 
projection of the NS orbit's semimajor axis, giving the masc 
function f(Mxf M^, i) = 4ir^(ax sin i)^/G P^orb = i/(l+q)2, 

where q = Mj^/My, i = angle between the planes of the sky and the 
orbit. ii) yields ay sin i of the visual orbit, and q = sin i/ 
Sjj sin i = Ky Porb (1-e )'2 /2ir a^ sin i, Ky = semiamplitude of the 
visual doppler velocity curve, e = eccentricity. iii) yields Sg, 
the eclipse half angle, linked to the stellar and Roche radii R* 
and Rl by R* = D (cos^ i + sin^ i sin^ 6e)^^2 = gD(a-b log q), where 
D = separation, a,b = constants depending on the ratio 

= “spin/“orb» ® ” fraction of Rj^ filled out by R*. Theoretical 
guesses at 6 and n lead then to a determination of Mj^, i and e. 

Both the observational errors in aj^ sin i, and 6g, and the 
guesses of ft, 6 lead to fairly large error bars. Values of are 
plotted in Figure 4, where results for 5 X-ray pulsars, as well as 
that of the binary radio pulsar are plotted. For more details, see 
(36,5). The remarkable thing is that, within these lo error bars, 
they all agree v;ith the range of masses (hatched area) expected on 
the basis of conventional nuclear and particle many-body physics. 

Related to mass determinations are also a number of obser- 
vations of a line feature at 400-460 keV in y-ray burst sources 
(37,38,39), interpreted as 511 keV e+e“ annihilation, redshifted by 
the gravitational field of a neutron star. Whi]^ the energy source 
need not be accretion, the redshift values provide information on 
the ratio M/R which are in agreement with expected values for 
neutron stars. This is true also of the March 5, 1979 burst (12), 
which showed pulsations indicative of being a neutron star (13). 
Calculations of annihilation spectra (40,41) indicate that the 
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"line" frequency shifts to the blue with increasing temperature, 
above a few x 10^ K, so the interpretation of the line frequency 
is not just given by gravitational redshift but the physics of the 
line formation as well. 

VI. MAGNETIC FIELD STRENGTH 

This is the area where most of the excitement has been lately. 

While initially researchers had to rely upon flux-freezing 
arguments, and use polarization measurements and theoretical 
considerations in radio pulsars to conclude that B ~ IO^^q^ a 
ground breaking direct experimental determination of the magnetic 
field strength of the X-ray pulsar Her X-1 came forth in 1977 
(42,43). A high altitude balloon experiment with a scintillation 
detector yielded a well defined line around 54 keV, with 15 o 
signi.ficance, which could only be interpreted as a cyclotron line 
in a field of B = 4.4 10^2 or 3.5 10^2 q if absorption (see 
Figure 3). This line was seen in three different flights by the 
original group, in one of the flights some evidence for a second 
harmonic being present (43), which however subsequently was not 
seen again, although this is consistent with the changes of the 
overall continuum between observations. While consistent results 
from other groups using similar resolution (~ 25 percent) (44,45) 
have been reported, a high resolution ( ~ 3 percent) Germanium 
experiment (46) came up with a feature at either 44 keV (emission) 
or 35 keV (absorption), and a line and continuum intensity which 
was within a factor 2 of previous values. Clearly further high 
resolution and high sensitivity measurements would be desirable to 
establish whether the line energy and intensity is a function of 
time, and perhaps of phase (e.g., (45)). Another X-ray pulsar for 
which a line feature has been reported is 4U0115+63 (47), and more 
recently, about thirty y-ray burst sources have been found to show 
well defined line features in the range 40-60 keV (39). implied 
field values would be in the range 2 10^2 1 b ~ 5 10^2 q, while 
theoretically there is no difficulty in understanding qualitatively 
the production of cyclotron line features in steadily accreting 
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X-ray pulsar atmosplieres (see below), the situation is less clear 
in the burst sources. Where relativistic temperatures exist 

at least initially, whereas the line widths suggest a lower 
temperature . 

VII. ‘ PHYSICS OP THE STRONGLY MAGNETIZED ATMOSPHERE 

On the theoretical side, this is the field which is most active at 
present. While nonmagnetic accreting neutron star atmospheres had 
been considered in (48,49,50), early work on magnetized atmospheres 
(51,52) had pointed out several interesting possibilities, such as 
cyclotron line features and intrinsic beaming. It was not however 
until the experimental detection of a cyclotron feature in Her X-1 
(42,43) that work on the subject began on a wide scale. Early 
attempts at interpretation (53,54,55,56) drove home the fact that 
progress could only be made after a number of basic questions on 
the elementary proceses were cleared up. 

The radiative cross sections were known to be strongly angle and 
frequency dependent in a magnetic field (57,58). Useful expres- 
sions were worked out in (60,61,62), elucidating the role of the 
ellipticity parameter and the normal polarization modes. 

Meanwhile, it had been pointed out that the vacuum polarization due 
to the virtual pairs in a strong magnetic field would influence 
polarization measurements (62). It was then realized that the 
vacuum polarization would also dominate the radiative cross 
sections and influence the spectrum, because of their dependence on 
the polarization normal modes (63,64,65). This effect can be 
thought of as being due to coherent forward scattering by virtual 
pairs, v^ich for real densities n^ ~ 10^"^ cm"^ and B ~ 10^^ G is 
more important than scattering off the real electrons. While in a 
cold e“ plasma, which rotates right handed if looked along only 
the right handed normal mode resonates at the cyclotron frequency, 
in the virtual e‘*'e“ plasma the right and left modes both 
resonate. This has also the effect of driving the usually elliptic 
modes towards being quasi linear. The large effect of the vacuum 
on the refractive index, and therefore on the polarization eigen- 
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inodes e ^((i>,6) of the plasma, has an equally large Impact on all 
the radiative cross sections, because these are expressions of the 
form do s <e |0|e>. The angle and frequency dependence of the 
Compton, bremsstrahlung and cyclotron processes are strongly 
affected. Also, a second resonance feature appears at the 
frequency » 3(4 x lO^^G/B) (n^/lO^^cm"^ heV, wher6 the real and 
virtual particles cancel their effects mutually (66,67). It also 
became apparent that "real” cyclotron absorption, where a photon 
disappears, its energy going into thermal motion, would not be 
important since radiative deexcitations would far exceed 
collisional ones at the resonance. The dominant process at is 
just resonant scattering (29,30). Collision of electrons with 
thermal protons (bremsstrahlung) would be the photon production 
mechanism, expressions being given in (29), based on (68). 


Using these cross sections, radiative transfer calculations were 
performed for uniform atmospheres, in the diffusion and two-stream 
approximations (29,30,66), including the frequency and polarization 
dependence of the cold plasma cross sections. These calculations 
revealed the importance of the mode exchange scattering which 

plays a crucial role in producing the cyclotron line (which would 
be absent if neglected). Photons are preferentially 

created in one mode, and escape in the other, depending on 
frequency. Incoherent scattering neglecting polarization effects 
was treated in (34) and in Monte Carlo calculations which also 
included thermal effects (31,32). Under coherent scattering, the 
cyclotron feature can appear either in emission or absorption, 
depending on whether a cooler spectrum is shined through a hotter 
gas or vice-versa (30). The inclusion of thermal and vacuum 
effects (69) can produce emission line profiles not unlike those 
observed in Her X-1, but this does not settle the emission versus 
absorption issue, until a self-consistent temperature structure for 
the atmosphere is known. The inclusion of quantum recoil effects 
in the cross sections and an‘ analysis of the redistribution 
function for one (resonant) polarization (70) raised the 
possibility of the line being asymmetrical respect to the rest 
resonance frequency. Thermal and quantum recoil effects were self- 
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consistently included in (35), where the effect of incoherent 
scattering inside and outside the resonance was investigated by 
means of a two-angle multi frequency Feautrier scheme. Incoherent 
scattering may explain the continuation of a power law continuum 
down to the ~ 1 keV range, but the photon input source function is 
a function of the still unsolved temperature and density structure 
of the atmosphere. 

The Coulomb collision of protons and electrons in a strong field 

had been investigated by (71,72,73) and, aside from possible 

collective effects (e.g., shock heating) this is the main candidate 

for a heating agent. Also, especially at low luminosities, where 

radiation pressure is negligible, and if a collionless shock does 

not arise, the stopping length of the infalling protons is a 

measure of the effective depth of the emitting atmosphere. The 

proton-electron binary collision cross section has been reevaluated 

in (74), where the effect,, of ion sound turbulence in the background 

plasma was considered, in (75), where knock-on collisions were also 

investigated, and in (76), where the e-e cross section was also 

evaluated. Except for (74), these results indicate longer stopping 

lengths than in the nonmagnetic case by a factor l/sin^ e, where 6 

is proton pitch angle, due to the fact that the electrons can only 

take up longitudinal momentum. Incoming fast protons at 6 ~ 0° 

pitch angle would in this case be stopped rather by nuclear p-p 

collisions, a ~ 50 g/cm . On the other hand, in (77,78) it is 

argued that, because the proton gyroradius is smaller than the 

Debye length of the plasma, the Coulomb stopping is more efficient 

2 

than previously thought, values of - 5-20 g/cm being put 
forward. The cross sections (75,76) were also used in investi- 
gating the electron distribution function (79), finding some 
evidence for departure from a Maxwellian. These studies of the 
electron heating rate, together with previously mentioned ones on 
the radiative transfer, should be combined to obtain a self- 
consistent structure for the X-ray atmosphere, necessary both for 
carculating a realistic spectrum, and a realistic beaming pattern. 
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Intrinsic beaming patterns have until noiy been calculated for 
guessed atmosphere parameters. Monte Carlo calculations for an 
optically thin cyclinder illuminated from below were done in (80) 
and (33 )« Whi!?h suggest a harder spectrum towards the equator. A 
plane parallel semi-infinite pure scattering atmosphere, with B 
perpendicular to the surface, was solved approximately in (81). 

This geometry leads to a pencil beam and can produce multiple 
pulses. The effect of absorption as well as scattering and 
polarization on the directionality was considered in 
82,83). In (19) the self-emission of a cylinder and a slab was 
solved by a Feautrier differential scheme, while in (18) the slab 
with background illumination and semi- infinite media with B 
perpendicular to the surface were solved with an integral scheme. 
Columns (fan beams) produce wider peaks than slabs (pencil beams) 
(19), and double peaks with a central hollow appear in pencil beams 
(18,19) as well as triple peaks and a hardening of the spectrum 
towards midpulse (18). These pulses appear to reproduce most of 
the major features of observed ones except the asymmetry and steep 
edges. To achieve this, it may be necessary to extend these 
calculations to distorted or non-dipole field configurations, 
although additional modulation by the Alfven surface sould also 
help. An intrinsic component seems however necessary to explain 
the multiple peak structure, the energy dependence of the pulse 
width and the phase dependence of the spectrum. 

The structure of the emitting atmosphere is largely still 
unknown. In the high luminosity pulsars, where the radiation 
pressure is sufficient to create a radiation-mediated shock, or 
deceleration zone, that stands off some distance dg, a "pill-box'' 
shaped atmosphere would arise at the base of the accretion 
column. Radiation would escape sideways (fan-beam) as well as 
perhaps upwards (85). The intriguing possibility also exists that, 
even in low luminosity pulsars, collisionless shocks may create a 
similar structure, but the plasma physics of such shocks is as yet 
unexplored in this context. Inhomogeneities may well arise due to 
instabilities (86). If shocks do not arise, the low luminosity 
pulsar atmospheres probably do not stick out significantly, the 



12 


1 2 • 
energetic (j nipV^^) proton b^am being decelerated by particle 

collisions in the denser part of the atmosphere. In this case, 

radiation escapes mostly upwards, producing a pencil bezun (52). A 

self-consistent atmosphere structure calculation has been done 

approximately only for the shock models (85,87). It would be very 

important to verify whether shocks, radiative or otherwise, indeed 

occur, and some simple tests may be possible. If the sides do not 

emit at the blackbody limit (<u> < wjj) the shock height dg ~ v tc, 

where V is post shock velocity and t^ is the Coulomb time for 

electron heating by the thermalized protons, so dg ~ n~^ ~ 

M ~ L . If the cyclotron line arises near the shock, one would 

expect a change of Ujj ~ (RnS + ds)”^ (i + const x L"i)“3. 

Similarly, the pulses, if produced by the fan beam, would become 

broader as L goes down, because the surface curvature sampled 

increases with dg, ~ ~ ~ shock 

arises, and one has an almost level hot spot producing a pencil 

beam, no such dependence of or on L would be expected, but a 

correlation should exist between the pulse multiplicity and the 

quantity For a wide range of i, a rough rule of thumb (18) 

is that pulses are single if ~ 10, triple if 10 ~ “h/*** ~ 

and double if ~ 50. One should point out that Her X-1 and 

4U0115+63, where u is known, have L > few 10^7 erg s-1, which 

probably leads to radiative shocks from which mostly fan beams are 

expected. Observations of such correlations or their absence could 

help settle the fan-pencil and shock-Coulomb controversy. 
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FIGURE CAPTIONS 

Figure 1. Relation between period decrement P and PL^/^, p being 
period, L the X-ray luminosity, for 9 X-ray pulsars, 
superimposed on theoretical curves based on disk 
accretion models (from Ref (7)). 

Figure 2. Pulse shapes and periods for a number of X-ray pulsars 
(from Ref. (5)). 

Figure 3. The X-ray spectrum of Her X-1, including the feature 
around 50 keV interpreted as due to the cyclotron 
effect in a field of 4.4 x 10 ^^ Gauss (from Ref. (43)). 

Figure 4. Neutron star mass limits derived from data on 5 X-ray 
pulsars and the binary radio pulsar. The shaded area 
indicates the range of values predicted by conventional 
equations of state (from Ref, (5)). 
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